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Abstract: Results from our recent investigations of quark number susceptibilities in 
both quenched and 2-flavour QCD are presented as a function of valence quark mass 
and temperature. A strong reduction (~40%) is seen in the strange quark susceptibility 
above in both the cases. A comparison of our isospin susceptibility results with the 
corresponding weak coupling expansion reveals once again the non-perturbative nature of 
the plasma up to 3Tc. Evidence relating the susceptibility to another non-perturbative 
phenomena, pionic screening lengths, is presented. 


1 Introduction 

Quantum chromodynamics (QCD) is the theory of interactions of quarks and gluons which 
are the basic constituents of strongly interacting particles such as protons, neutrons and 
pions. A complete lack of any experimental observation of a free quark or gluon led to 
the hypothesis of their permanent conhnement in the observable particles. As pointed 
out by Satz in his talk, relativistic heavy ion collision experiments at BNL, New York, 
and CERN, Geneva offer the possibility of counting the number of degrees of freedom of 
strongly interacting matter at high temperatures, thereby providing a strong argument for 
Heisenberg to accept the physical reality of quarks and gluons in spite of their confinement 
at lower temperatures. The basic idea here is to look for the production of quark-gluon 
plasma (QGP), predicted by QGD to exist beyond a transition temperature Tc, in the 
heavy ion collisions and when found, determine its thermodynamical properties, such as, 
e.g., its energy density. At very high temperatures, these properties can be computed 
theoretically and are seen to be directly proportional to the number of quarks and gluons. 
In order to test this idea, however, one has to face the fact that the temperatures reached 
in the current, and near future, experiments are likely to lie below about 5-10 Tc. The 
current best theoretical estimates for thermodynamic observables in this temperature 
range are provided by numerical simulations of QGD defined on a discrete space-time 
lattice 0. It seems therefore prudent to evaluate as many independent observables as 
possible and compare them with both experiments and approximate analytical methods. 
Quark number susceptibilities constitute a useful independent set of observables for testing 
this basic idea of counting the degrees of freedom of strongly interacting matter. 

^Based on work done ]ll, u\ with Sourendu Gupta and Pushan Majumdar. 




Investigations of quark number susceptibilities from first principles can have direct 
experimental consequences as well since quark flavours such as electric charge, strangeness 
or baryon number can provide diagnostic tools for the production of flavourless quark- 
gluon plasma in the central region of heavy ion collisions. It has been pointed out recently 
i,i on the basis of simple models for the hadronic and QGP phases that the fluctuations 
of such conserved charges can be very different in these two phases and thus can act as 
probes of quark deconhnement. Indeed, excess strangeness production has been suggested 
as a signal of quark-gluon plasma almost two decades ago [0. Lattice QCD can provide a 
very reliable and robust estimate for these quantities in both the phases since in thermal 
equilibrium they are related to corresponding susceptibilities by the fluctuation-dissipation 
theorem: 

{SQ^) oc = Xq(T ,fio = 0) . (1) 

Here fiQ is the chemical potential for a conserved charge Q, and Z is the partition function 
of strongly interacting matter in volume V at temperature T. Unfortunately, the fermion 
determinant in QCD becomes complex for any nonzero chemical potential for most quan¬ 
tum numbers including those mentioned above. Consequently, Lattice QCD is unable 
to handle hnite chemical potential satisfactorily at present, and cannot thus yield any 
reliable estimates of any number density. However, the susceptibility above, i.e., the hrst 
derivative of the number density at zero chemical potential, can be obtained reasonably 
well using conventional simulation techniques, facilitating thereby a nontrivial extension 
of our theoretical knowledge in the nonzero chemical potential direction. 

Quark number susceptibilities also constitute an independent set of observables to 
probe whether quark-gluon plasma is weakly coupled in the temperature regime accessible 
to the current and future planned heavy ion experiments (say, 1 < T /T^ < 10). A lot of the 
phenomenological analysis of the heavy ion collisions data is usually carried out assuming 
a weakly interacting plasma although many lattice QCD results suggest otherwise. It has 
been suggested [Q that resummations of the finite temperature perturbation theory may 
provide a bridge between phenomenology and the lattice QCD by explaining the lattice 
results starting from a few T^. As we will see below, quark number susceptibilities can 
act as a cross-check of the various resummation schemes. Earlier work on susceptibilities 
0 did not attempt to address this issue and were mostly restricted to temperatures very 
close to Tc- Furthermore, the quark mass was chosen there to vary with temperature 
linearly. We improve upon them by holding quark mass fixed in physical units [m/Tc = 
constant). We also cover a larger range of temperature up to 3 Tc and the accepted range 
of strange quark mass in our simulations. 

2 Formalism 

After integrating the quarks out, the partition function Z for QCD at finite temperature 
and density is given by 

Z = y det M(mu,/i„) det /id) det M(m 5 , ps). (2) 

Here {Uy{x)}, v = 0-3, denote the gauge variables and Sg is the gluon action, taken to 
be the standard 0 Wilson action in our simulations. Due to the well-known “fermion 
doubling problem” [Q, one has to face the choice of fermion action with either exact 
chiral symmetry or violations of flavour symmetry. The Dirac matrices M depend on this 



choice. Since we employ staggered fermions, the matrices, M, are of dimensions SN^Nt, 
with Ns{Nt) denoting the number of lattice sites in spatial(temporal) direction. These 
fermions preserve some chiral symmetry at the expense of flavour violation. Although, 
they are strictly dehned for four flavours, a prescription exists to employ them for arbitrary 
number of flavours which we shall use. mj and fif are quark mass and chemical potential 
(both in lattice units) for flavour /, denoting up(u), down(d), and strange(s) in eq.(|D. 
The chemical potential needs to be introduced on lattice as a function g{fi) and g{—g) 
multiplying the gauge variables in the positive and negative time directions respectively, 
such that [|T^] i) g{fi) ■ g{—g) = 1 and ii) the correct continuum limit is ensured. While 
many such functions g can be constructed, exp(/i) being a popular choice, the results for 
susceptibilities at /x = 0 can easily be shown to be independent of the choice of g even 
for hnite lattice spacing a. From the Z in eq. (|]), the quark number densities and the 
corresponding susceptibilities are dehned as 
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where the volume V = and the temperature T = To lighten the notation, 

we shall put only one subscript on the diagonal parts of x- 

In order to obtain information for quark-gluon plasma in the central region, we evaluate 
the susceptibilities at the point fif = 0 for all /. In this case, each Uf vanishes, a fact that 
we utilize as a check on our numerical evaluation. Moreover, the product of the single 
derivative terms in eq. (|^) vanishes, since each is proportional to a number density. We 
set TTiu = rrid < rris- Noting that staggered quarks have four havours by default, Nf = 4, 
and dehning fi^, = fiu — one hnds from eq. that the isotriplet and strangeness 
susceptibilities are given by 


X3 = Xs = ^[Oi{ms) + , (4) 

where C>i = (Tr(M"M-i - O 2 = {{TiM'M' = dM/dfi and 

M" = d^M/dfi^. The angular brackets denote averaging with respect to the Z in eq. (^. 
One can similarly dehne baryon number and charge susceptibilities. We refer the reader 
for more details on them to Ref. [^. 

In the discussion above, quark mass appears as an argument of Oi and implicitly in 
the Boltzmann factor of Z. Let us denote it by m^ai and nisea respectively. While the 
two should ideally be equal, we evaluated the expressions above in steps of improving 
approximations (and increasing computer costs) by hrst setting m^ea = cc for all flavours 
(quenched approximation 1 ) and then simulating two light dynamical flavours, by setting 
nisea/Tc = 0.1 (2-flavour QCD Q). In each case we varied m^ai over a wide range to cover 
both light u,d quarks as well as the heavier strange quark. Details of our simulations 
as well as the technical information on how the thermal expectation values of Oi were 
evaluated are in Refs. |l|, . 


3 Results 

Based on our tests |ll] of volume dependence, made by varying W from 8 to 16, we chose 
Ns = 12, as the susceptibilities differed very little from those obtained on an W = 16 
lattice. Choosing Nt = 4, the temperature is T = (4a)where the lattice spacing a 
depends on the gauge coupling f3. Using the known PciNj.) for Nj. = 6, 8,12, where j3c is 












Figure 1: xIxfft as a function of temperature for various valence quark masses. 


the gauge coupling at which chiral (deconhnement) transition/cross-over takes place for 
lattices with temporal extent Af/, we obtained results at T/Tc = N'^/Nt =1.5, 2 and 3 in 


both the quenched approximation and the 2-flavour QCD. From the existing estimates 
of Tc for 2-flavour QCD, one hnds that the sea quark mass in our dynamical simulations 
corresponds to 14-17 MeV. Fig. ^ shows our results, normalized to the free held values on 
the same size lattice, Xfft- These can be computed by setting gauge helds on all links 
to unity: Uu{x) = 1 for all v and x. Note that due this choice of our normalization, the 
overall factor Nf for degenerate Favours cancels out, permitting us to exhibit both X 3 and 
Xs of eq.(§) on the same scale in Fig. The continuous lines in Fig. were obtained from 
the interpolation of the data obtained in the quenched approximation (the data are not 
shown for the sake of clarity), while the results for two light dynamical Favours are shown 
by the data points. In each case the value of rriyai/Tc is indicated on the left. Due to the 
fact [0 that the contribution of O 2 to Xs turns out to be negligibly small for T > Tc, X3 
and Xs appear coincident in Fig. For the real world QCD, the low valence quark mass 
results are relevant for xs and those for moderate valence quark masses are for Xs- 

Although Tc diFers in the quenched and 2-Favour QCD by a factor of 1.6-1.7, the 
respective susceptibilities shown in Fig. as a function of the dimensionless variable 
T/Tc change by at most 5-10% for any rrivai/T. Thus the eFect of “unquenching”, i.e., 
making 2 Favours of quarks (u and d) light enough to include the contribution of the 
corresponding quark loops, appears to be primarily a change of scale set by Tc- Since the 
strange quark is a lot heavier than the up and down quarks, this suggests further that 
including its loop contributions, i.e., including a dynamical but heavier strange quark, may 
not change the results in Fig. |I| signiFcantly. For a wide range for strange quark mass 
of 75 to 170 MeV, the strangeness susceptibility can be read oF from the shaded region. 
It is smaller by about 40% compared to its ideal gas value near Tc and the suppression 
shows a strong temperature dependence. This has implications for the phenomenology of 
particle abundances, where an ideal gas model without such a suppression of strangeness 
is employed and could therefore result in underestimates of the temperatures reached. 

For the smallest rriyai/T and highest temperature we studied, the ratio x/Xfft in 
Fig. g is seen to be 0.88 (0.85) for 2-Favour (quenched) QCD, with a mild temperature 
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Figure 2: xIxfft as a function of Os/tt for Nf dynamical massless quarks. 


variation in the large T-region. Since its variation with valence quark mass is negligibly 
small for small rrivai/T, one can assume the results for massless valence quarks to be 
essentially the same as those for m„ai/T = 0.03 in Fig. |l|. In order to check whether the 
degrees of freedom of QGP can really be counted using these susceptibilities, one needs 
to know whether the deviation from unity can be explained in ordinary perturbation 
theory or its improved/resummed versions. Usual weak coupling expansion 0 yields 

x/Xfft = 1 — 2^[1 + ^)]. Fig. 1^ shows these predictions for various Nf along 

with the leading order A^j-independent prediction. Using a scale 2ttT for the running 
coupling and Tc/Ajjg = 0.49(1.15) for the Nf = 2(0) theory |]^, the values T/Tc =1.5 
and 3 are marked on the hgure as the second (hrst) set. As one can read off from the 
Fig. 0, the ratio decreases with temperature in both the cases in the range up to 3Tc 
whereas our results in Fig. |l] display an increase. Furthermore, the perturbative results 
lie signihcantly above in each case, being in the range 1.027-1.08 for 2-flavour QCD and 



Figure 3: (open symbols) and x,r/10T^ (filled symbols) as a function of M^^/T at 2Tc (circles) 

and 3Tc (boxes). 












0.986-0.994 for quenched QCD. Although the order of magnitude of the degrees of freedom 
can be gauged from these results and their eventual comparison with experiments, they 
do call for clever resummations of perturbation theory for a more convincing and precise 
count. 

Alternatively, the deviations from free field theory could stem from non-perturbative 
physics. One known indicator of non-perturbative physics in the plasma phase is the 
screening length in the channel with quantum numbers of pion. While it exhibits chiral 
symmetry restoration above by being degenerate with the corresponding scalar screen¬ 
ing length, its value is known to be much smaller than the free field value unlike that for 
other screening lengths. Fig. || shows our results for ys and Xtt (defined as a sum of the 
pion correlator over the entire lattice) as a function of the inverse pionic screening length, 
Mt,/T. It suggests the non-perturbative physics in the two cases to be closely related, if 
not identical. 
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